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In the present work we study the well-known Two Capacitor Problem from a new perspec- 
tive. Although this problem has been thoroughly investigated, as far as we know there are 
no studies of the thermodynamic aspects of the discharge process. We use the Free Electron 
Gas Model to describe the electrons' energy levels in both capacitors in the low temperature 
regime. We assume that the capacitors and the resistor can exchange energy freely with a 
heat reservoir. We assume that the resistance is large enough to consider an isothermal heat 
exchange between the resistor and the heat reservoir. Thereby we obtain a positive entropy 
variation due to the discharge process, corroborating its irreversibility. 
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I. INTRODUCTION 

In order to compute the entropy change AS, suppose that the circuit of Fig. (OQ) is in thermal 
contact with a heat reservoir at temperature Tq. Consider a simple circuit that consists of two 
identical capacitors, with the same capacitive constant C, and one resistor R, all connected in 
series. One of the capacitors (capacitor 1) is initially charged with a charge qo, while the other 
one has zero charge (capacitor 2). The circuit also has a switch that prevents the flow of current, 
as shown in figure (OQ). 
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FIG. 1. Illustration of the circuit described in the text. In (a) the switch is open and there is no current 
flowing. In (b), we have the final equilibrium situation, in which a sufficiently large time has elapsed for 
the system to reach equilibrium. 



If we close the switch a current will flow until the two capacitors reach the same potential 
V (in our case, V = since the capacitors have the same capacitive constant C). We know 
intuitively, that the discharging process is highly irreversible: the system reaches spontaneously 
the equilibrium configuration described above and it will not go back to the initial configuration 
without external interference. There is a simple analogy between this discharging process and 
the isothermal expansion of an ideal gas (commonly used to obtain the entropy change in a free 
expansion). In both cases, there is a positive entropy change for the gas, due to the redistribution 
of particles. However, here the particles in the play are fermions, thus obeying a different 
statistical distribution, which is a fundamental difference between then. 

There are many previous works which deal with this simple system using a variety of ap- 
proaches. In iu,\2\, the authors focus on energy considerations, while |3|-l50 discuss specifically 
the electromagnetic radiation produced by the discharge process. Some of those works neglect 
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the circuit's electrical resistance, but include a self inductance. Some authors consider both 
resistance and inductance. 

In [|6|] the author measures some quantities in order to obtain the entropy change of the 
capacitor charging process. The author shows experimentally that if the charging process can 
be considered a quasistatic one, the entropy change due to resistor and the heat reservoir adds 
to zero. In our discharging problem, this limiting case is achieved by a a large resistance R. 
However, this work considers only the entropy variation due to heat exchange between the 
resistance and the heat reservoir, obtaining AS = when the process is sufficiently slow. We 
will use a simple model to discuss the entropy variation of the whole system. 

It is our goal to calculate the change in the thermodynamical entropy S due to the discharge 
process. We will show how this can be done using a simple model to describe the electrons in the 
capacitors, allowing a straightforward evaluation of AS. As mentioned before, the discharging 
process is clearly irreversible, such that any model chosen to describe the discharging process 
must provide AS > 0. 

H. ENTROPY CALCULATION 

According to the second law of thermodynamics, any process that occurs in an isolated 
system must have AS > 0. So, to compute AS, we need the entropy change of all components 
of our system: both capacitors, the heat reservoir and the resistance. First, we will evaluate AS 
for the heat reservoir and the resistance. 

As mentioned above, we consider that the resistance is large enough, in order that there is a 
isothermal heat exchange between the resistor and the heat reservoir. Thus AS^ = — AShr and 
the total entropy change due to this exchange is zero. As a consequence the total entropy varia- 
tion of the discharge process is only due to the charge rearrangement in the two capacitors. To 
calculate this variation we will use the Free Electron Gas model to describe the electrons in the 
metal plates of the capacitors. We can consider the plates as three dimensional: although their 
thicknesses are very small in comparison to the other dimensions, they are indeed macroscopic 
as far as the motions of electrons are concerned. 

The Free Electron Gas Model was the first attempt at a microscopic description of the prop- 
erties of metals. It was first used by Drude[7] and later modified by Sommerfeld O, |8|] to 
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include the quantum nature of the electrons. The model consists of the assumption that elec- 
trons are confined to a box with impenetrable walls. Inside the box the electrons are subject 
to a uniform potential Vq. Although this model is very simple, it describes qualitatively very 
important features of metals, such as the Wiedemann-Franz Law and the fact that the electronic 
contribution to the specific heat is proportional to T at low temperatures (as compared to the 
Fermi temperature) QV|, LSQ . 

Since the model assumes that electrons do not interact with each other, the solution to the 
problem consists of finding the eigenstates of a quantum particle confined to a box and filling 
up these states according to Fermi-Dirac statistics. In order to obtain the expression for the 
electronic contribution to the entropy, we will use a well-known result for the specific heat in 
this model[ 7, <!), 

Nk B n 2 k B T 

where k B is the Boltzmann constant and £f is the Fermi energy which is given by 

e -^i-^(^L) 2/3 N 2/3 (2) 
2m 2m \ V ) 

where N is the total number of free electrons, h is the Planck constant divided by 2%, m is 
the electron mass and V is the volume of the macroscopic sample. The expression (OQ) results 
from the Sommerfeld expansion for the specific heat, which is very accurate as long as the 
temperature is small when compared to the Fermi temperature Tf. This is easily satisfied at 
room temperature, at which T/Tp ~ 10 ). 
Using the following thermodynamic identity, 

Cv = t(^) , (3) 



dT 



v 



it is clear that the entropy S is equal to the specific heat Cy for the Free Electron model II 1211 . 

^_ Nk B K 2 k B T 

Since we assume that the discharging process occurs at constant temperature, we can rewrite 
the Eq. © as 

S(T,V,N)=A(T,V)N l /\ (5) 
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where 

2/3 

(6) 



ACr.v)^^-^ r>o. 



Now that we have the entropy as a function of all the relevant parameters of the problem, 
we are able to evaluate the entropy change of the system due to the discharge process. We 
denote the entropy of each capacitor i in the initial / or final F configuration by by Sj and 
Sp) respectively. The total entropy change is given by the sum of the entropy changes of each 



capacitor, 



AS = ASi +AS 2 = (S { p ] -Sf } ) + (4 2) -Sf } ). (7) 



Before we go on performing the calculations, let's make an observation about the Free Elec- 
tron Gas model. We assume that the piece of metal that we are modeling possess N atoms, 
and that each one of these atoms contributes with a electrons to our sample (typically, a is 1 or 
2). When an atom donates a electrons, it will be positively charged, with q = ae, where e is 
the electron's fundamental charge. So, when we consider all the sample, we conclude that our 
piece of metal is uncharged because the charge of the electrons cancels out the charge of the 
ions, although there exists free electrons. 

Furthermore, we are considering metal plates that could be positively or negatively charged. 
We know that the excess charge q that a macroscopic metal can support is very small, when 
compared to the number of valence electrons in the sample. To see a good discussion on this 
subject, see [10], when Feynman introduces electric forces. 

With these observations in mind, let us continue with our computations. Consider first the 
two plates of capacitor 1. Initially, one plate has an excess charge qo, and the other, —qo. We 
known that the electric charge is quantized, so we have go = N^e, where Nq is the excess number 
of electrons in the plate. Remembering that the total number of electrons in a plate is Na plus 
the excess charge, we have 



ASi=A{T,V) 



Na + |) 1/3 + (Na - g) ^ - (Na + ^ - ( N a - f ^ 



(8) 



because in the equilibrium situation the excess charge in each plate is ±qo/2, as discussed 
before. 

Analogously, for the capacitor 2, we have 
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AS 2 =A(r,V) 



1/3 



(9) 



Now, to obtain AS, we substitute Eqs. ® and © into ©: 



AS = A(r,V)(Afa) 1/3 2(1+x) 1/j +2(1-jc 



,1/3 



,1/3 



(1 +2jc) 1/3 - (1 -2jc) 



1/3 



(10) 



where we defined x = qo/lNae, the ratio between the excess charge and the number of valence 
electrons in the neutral sample. The expression obtained for AS in (PTOl) is always positive, as 
shown in Fig. ©. 

Although we obtained a expression for AS, and have shown that AS > 0, we can go a little 
further. As we discussed earlier, typically we have x <C 1 . Thereby, we can expand the Eq. (flOl) . 
using the binomial expansion for all terms. Doing this, and collecting only second order terms, 
we obtained 



AS 



AA{T,V)x 2 (Na) 1 / 3 A{T,V)qftNa) 



-5/3 



(ID 



9 9e 2 
In this case, we needed to consider second order terms in the expansion for AS, because the first 
order term vanishes. This can be noticed visually in Fig.©, as the slope of the curve vanishes 
as x — > 0. 



III. ENERGY CALCULATION 



Up to now we have discussed the entropy change in the discharging process. We can also 
discuss what happens with the energy. First, we note a somewhat counter-intuitive fact. The 
electrostatic energy change AE E ; of the system, between the initial and the final configurations 
does not depend on the resistance R. As long as R > lllin . we have 



AE E i = - 



4o 
AC 



(12) 



However note that the eq.([T21) only takes into account the electrostatic interaction. In order to 
obtain the total energy variation we should consider the energy change due to the rearrangement 
of the energy levels of the free electron gas. At room temperature the total energy of the free 
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FIG. 2. Entropy change as a function of the ratio between the excess charge and the total number of 
valence electrons. 

electron gas is well approximated by its ground state energy given by 

E = ^Ne F . (13) 

Although in the entropic calculation we used the first order Sommerfeld expansion, in the 
energy calculation the use of the ground state energy is justified by the fact that in the latter 
AE ^ even for T = K. For the entropy, however, AS = for T = K and going beyond 
zeroth order is essential. 

Using eq. (fT3l) to calculate the energy variation of the four capacitor plates, we obtain [see 
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fig.©], analogously to the entropy calculation, 



AE FE = B{V)N 5/3 [2(1 +jc) 5/3 + 2(1 -jc) 5/3 - (1 +2jc) 5/3 - (1 - 2x) 5/3 - 2 



(14) 



where 



3 fh 2 \ /3t^ 2/3 



(15) 




FIG. 3. Energy change as a function of the ratio between the excess charge and the total number of 
valence electrons. 



Expanding (|T4l) for x <C 1 we obtain 
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AE FE = -—B{V)N 5 I\ 2 . 



(16) 



Thus, the total energy variation is given by 



± E = -^--B{V)N 5 l 3 x 2 . 
4C 9 v ' 



(17) 
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IV. CONCLUSIONS 



In this work we introduced a new approach for the two-capacitor problem. Instead of looking 
at it from an electromagnetic point of view, we focus on its thermodynamic properties, mainly 
the entropy change due to the discharge process. Our goal was to show that this process is 
irreversible, which means that the entropy variation must be always positive (AS > 0). We 
assumed that the resistance is in thermal contact with a heat reservoir at the same temperature, 
which leads to ASr + AShr = 0. Thus, all the entropy variation is due to the change in the 
electrons' distributions in the plates of the capacitors. We assumed that the valence electrons in 
the metal plates can be modeled by the Free Electron Gas Model. Under these considerations, 
we obtained AS > and AE < for the whole process. 
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